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The  present  paper  investigates  the  enzymatic  hydrolysis  of  tuna  fish  oil  with  Candida  rugosa  lipase  (CRL)  in 
biphasic  solvent  system  for  the  production  of  free  fatty  acids  (FFAs).  Effect  of  reaction  parameters  were  studied  by  varying 
suitable  reaction  conditions  such  as  pH,  temperature,  agitation  speed,  water  and  solvent  concentrations.  When  the  reaction 
was  carried  out  for  24  h,  86.5%  hydrolysis  was  achieved.  For  hydrolysis  kinetics,  the  model  given  by  Prazeres  et  al.,  [1] 
with  second  product  inhibition  was  used.  The  kinetic  model  was  fitted  using  MATLAB®  to  determine  the  best  kinetic 
parameters.  R2  and  root  mean  square  error  (RMSE)  were  found  to  be  0.961  and  56.7,  respectively.  The  average  value  of 
kinetic  constants  using  the  Prazeres  model  were  estimated  as  KM=  4.26  (J.  moles  FFA/ml,  K;,  =  6.0  X  10"6  (J.  moles  FFA/  mg 
enzyme-  h,  Ki2  =  0.042  \x  moles  FFA/  mg  enzyme  h  and  K2  =  122.3  (i  moles  FFA/  mg  enzyme-  h. 
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The  importance  of  polyunsaturated  fatty  acids  (PUFA)  in  human  nutrition  and  disease  prevention  was 
scientifically  recognized  three  decades  ago.  In  today's  perspective,  the  social  demand  is  to  focus  on  the  synthesis  of 
compounds  which  can  improves  the  immunity  against  several  chronic  and  degenerative  diseases  like  dyslexia,  dyspraxia, 
alzheimer,  schizophrenia  and  heart  diseases  [2-4].  Omega-3  long  chain  PUFA  such  as  DHA  (docosahexaenoic  acid)  is 
widely  reported  to  have  pharmaceutical  and  neutraceutical  properties.  DHA  is  not  only  essential  for  the  development  of 
human  brain  but  also  used  in  the  prevention  and  control  of  various  human  diseases  and  disorders  such  as  cardiovascular 
disease,  inflammation,  allergy,  cancer,  immune  response,  diabetes,  hypertension  and  renal  disorders  [5].  Tuna  fish  oil  is 
one  of  the  richest  natural  sources  of  DHA  with  its  concentration  in  the  range  of  25-30  (wt  %).  The  extraction  of  DHA  from 
fish  oil  has  been  reported  with  both  chemical  and  enzymatic  hydrolysis  [6-8].  Enzymatic  hydrolysis  occurs  with  high 
specificity  at  ambient  temperature  and  pressure  in  comparison  with  the  conventional  fat  splitting  process  [9]  and  therefore 
preferred  over  conventional  methods  from  last  few  decades  [10,  11].  The  advantages  of  enzymatic  processes  over 
conventional  chemical  processes  include  mild  reaction  conditions  of  pH,  temperature,  reduced  energy  cost,  low  amount  of 
enzyme  consumption,  high  catalytic  efficiency  and  wide  range  of  fatty  acid  selectivity  of  lipases  which  leads  to  specific 
and  pure  products  [12,  13].  This  makes  the  enzymatic  process  much  cleaner  and  energy  efficient  than  the  conventional 
thermal  fat  splitting. 

Lipases  (commonly  known  as  triacylglycerol  acylhydrolyase  and  classified  as  EC  3.1.1.3)  are  used  industrially 
[14-17]  and  commercially  in  the  breakdown  of  the  acyl  glycerol  bond  in  the  presence  of  the  sufficient  quantity  of  the  water 
and  lead  to  the  formation  of  the  monoglycerides,  diglycerides,  free  fatty  acids  and  the  glycerol  as  by-products  [18,19]. 
Lipases  have  gained  the  attention  of  the  researchers  worldwide  because  of  their  applications  in  various  streams  such  as 
synthesis  of  oleo-chemicals,  fine  chemicals,  pharmaceuticals,  fragrances  and  flavor,  value  added  food  products  with 
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medicinal  values  etc.  [2,  3].  Lipases  are  interfacial  enzymes  remained  activated  at  the  oil-water  interface  due  to  the  open 
lid  conformation  of  lipases  at  the  interface  of  biphasic  system  [20,  21].  The  mechanism  involves  the  attachment  of  acyl 
group  (RCOO  )  on  the  positive  side  chain  (NH3+)  and  hydrogen  ion  (H+)  on  the  negative  side  chain  (COO  )  on  the  lipase 
active  site.  This  causes  the  release  of  free  fatty  acids  (RCOOH)  in  the  solvent  biphasic  system  and  lipase  can  also  be 
recovered  after  the  extraction  of  solvent  phase. 

The  kinetics  of  enzymatic  hydrolysis  using  lipases  has  been  studied  primarily  in  monolayer  by  Verger  et  al.  [22] 
assuming  that  all  the  products  of  reactions  were  soluble  in  the  water  phase,  diffuse  rapidly  away  from  the  interface  and 
induce  no  change  in  substrate  with  time.  The  above  assumptions  are  not  true  for  practical  cases  such  as  hydrolysis  of 
oils/fats  containing  long  chain  fatty  acids.  Martinez  et  al.  [23]  applied  the  kinetic  model  of  Verger  et  al.  [22]  for  hydrolysis 
of  tributyrin  by  Candida  cylindracea  lipase  to  derive  a  rate  equation  as  a  function  of  volume  fraction  of  dispersed  phase,  <t>, 
and  the  speed  of  agitation,  N.  They  used  the  equation  of  Tavlarides  and  Bapat  [24]  to  calculate  the  mean  drop  diameter  of 
dispersed  phase.  They  verified  the  model  as  adequate  and  found  the  influence  of  agitation  speed  low  and  limited  to  the 
domain  of  low  O  (<E>  <  10%).  Tsai  et  al.  [25]  studied  the  kinetics  of  enzymatic  hydrolysis  of  olive  oil  in  biphasic  organic- 
aqueous  systems  considering  the  reaction  as  single  substrate  reaction.  The  model  does  not  talk  about  substrate  and  product 
inhibition  and  hence  the  validity  of  this  model  in  long  term  hydrolysis  is  doubtful.  Hermansyah  et  al.  [26]  gave  Ping  Pong 
Bi  Bi  model  for  the  triolein  hydrolysis  using  Candida  rugosa  lipase  in  the  biphasic  oil-water  system.  The  model  described 
the  effect  of  the  initial  enzyme  concentration,  the  interfacial  area  and  the  initial  concentrations  of  triolein  and  a  fatty  acid 
on  the  entire  process  of  the  stepwise  triolein  hydrolysis. 

Brijwani  and  Vadlani  [27]  proposed  a  kinetic  model  for  hydrolysis  of  corn  distillers  dried  grain  with  solubles 
(DDGS)  oil  in  emulsion  system  which  uses  the  effect  of  interfacial  area  and  product  inhibition  terms,  both  simultaneously 
at  higher  agitation  speed  and  high  enzyme  concentration  and  have  also  reported  that  at  low  enzyme  concentration  and  low 
speed  of  agitation,  no  inhibition  was  observed. 

Knezevic  et  al.  [28]  reported  lipase  catalyzed  hydrolysis  of  palm  oil  in  a  lecithin/iso -octane  reverse  micellar 
system  followed  one  substrate  first-order  reversible  Michaelis-Menten  kinetics,  when  the  initial  substrate  concentration  is 
less  than  0.325  mol/dm3.  Yao  et  al.  [29]  studied  the  kinetics  of  lipase-catalyzed  hydrolysis  of  olive  oil  in  AOT/isooctane 
reversed  micellar  system  and  proposed  a  model  incorporating  product  inhibition,  enzyme  deactivation  and  influence  of 
both  substrate  and  product  on  enzyme  stability  in  reverse  micellar  system. 

Prazeres  et  al.  [1]  proposed  a  kinetic  model  with  nonlinear  third  order  product  inhibition  and  found  it  to  be  most 
adequate  for  the  description  of  long  term  reaction  kinetics  during  batch  hydrolysis  of  olive  oil  with  Chromobacterium 
viscosum  lipase  in  AOT/isooctane  reverse  micellar  system.  The  model  proposed  by  Prazeres  et  al.  [1]  was  single  substrate 
model  as  the  water  content  used  in  their  reaction  system  was  very  high.  It  can  be  seen  that  the  various  kinetic  models  given 
in  the  literature  are  based  on  assumption  of  one  substrate  reaction  and  some  of  them  also  incorporated  product  inhibition 
and  effect  of  other  reaction  conditions  like  interfacial  area  and  agitation  speed  etc.  However,  it  is  well  known  that  in  the 
hydrolysis  reaction,  two  substrates  take  part  in  reaction,  and  therefore  the  reaction  may  be  considered  as  second  order 
reaction.  When  large  excess  of  one  reactant  is  taken,  then  the  reaction  may  be  considered  as  pseudo  first  order  reaction. 

In  the  present  work,  the  fitting  of  experimental  data  with  the  kinetic  model  given  by  Prazeres  et  al.  [1], 
incorporating  second  order  product  inhibition  has  been  studied.  Also,  the  agitation  speed  during  the  reaction  has  been  kept 
sufficiently  high  where  the  effect  of  interfacial  area  at  the  given  concentration  of  enzyme  becomes  negligible  and  the  water 
concentration  has  been  kept  sufficiently  high. 
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MATERIALS  AND  METHODS 
Materials 

Tuna  fish  oil  of  refined  grade  was  purchased  from  Jedwards  International  Inc.,  MA,  U.S.A.  Lipase  enzyme  such 
as  Candida  rugosa  has  random  positional  specificity  with  4.01  U/mg  activity  was  supplied  by  Sigma  Aldrich,  Canada.  Iso- 
octane,  sodium  methoxide,  10%  solution  of  BF3  in  methanol,  DHA,  EPA,  fatty  acid  methyl  ester  standards  used  for  GC 
calibration  and  Triolein,  Diolein,  Monolein  used  for  HPLC  calibration  were  purchased  from  Sigma  Aldrich,  Ontario, 
Canada.  All  other  chemicals  like  potassium  hydroxide  (KOH),  ethanol,  hexane,  anhydrous  sodium  sulfate,  hydrogen 
chloride  (HO)  and  methanol  were  obtained  from  VWR,  Saskatoon,  Canada. 

Experimental  Procedure 

The  suitable  reaction  conditions  for  hydrolysis  of  tuna  fish  oil  with  CRL  was  studied  at  different  pH,  temperature, 
agitation  speed  (rpm),  oil  to  water  and  oil  to  solvent  ratio.  Experiments  were  carried  out  in  triplicates  at  different  reaction 
parameters  by  varying  only  one  reaction  parameter  at  a  time  and  keeping  other  parameters  constant.  The  reaction  mixture 
consisted  of  lg  fish  oil  with  1ml  of  isooctane,  40  mg  lipase  enzyme  dissolved  in  10  ml  distilled  water  was  kept  in  a  50  ml 
conical  flask  and  agitation  was  carried  out  using  a  magnetic  stirrer.  Methanol  was  added  to  stop  the  reaction  and  acid  value 
of  the  reaction  mixture  was  measured  by  titration  with  0.5N  KOH  to  determine  the  formation  of  FFAs  and  the  percentage 
hydrolysis  was  calculated  by  the  following  formula, 

Acid  Value  (hydrolysed  oil  —  Blank') 

%  Hydrolysis  =  —  —  — -xl00% 

Saponification  Value  (original  oil)  —  Acid  value  (original  oil) 

Standard  AOCS  test  method  was  used  for  estimating  the  acid  value  [30]  of  collected  samples  at  different  reaction 
time  to  determine  amount  of  free  fatty  acids  (FFAs)  formed.  The  activity  of  lipase  has  been  expressed  as  pmoles  of  FFAs 
formed  per  ml  of  reaction  mixture.  At  the  suitable  conditions,  reaction  was  carried  out  for  24  h  and  the  kinetic  study  was 
done  by  changing  the  substrate  concentrations. 

The  amount  of  oil  was  varied  in  the  range  of  0.8  g  to  1.5  g  (0.8  g,  1.0  g.  1.25  g,  1.5  g).  Maintaining  the  batch 
volume  constant  i.e.  12  ml,  40  mg  CRL  was  dissolved  in  10  ml  water  while  changing  the  iso-octane  amount  from  0.5  to  1.2 
ml.  For  the  kinetic  study,  individual  batch  of  experiments  have  been  carried  out  at  initial  substrate  (oil)  concentrations  of 
214  (Soi),  267  (S02),  334  (S03)  and  401  (S04)  |jmoles/ml,  for  different  reaction  time  and  FFAs  produced  have  been 
determined.  Further,  MATLAB®  software  was  used  to  determine  the  kinetic  parameters.  R2  and  Root  mean  square  error 
(RMSE)  were  evaluated  for  determining  the  goodness  of  fitting. 

Analysis:  The  fatty  acid  composition  of  tuna  fish  oil  was  determined  by  using  an  Agilent  make  Gas  Chromatography 
system  (model  5890A),  equipped  with  flame  ionization  detector  (FID,  250°C)  and  capillary  column  DB-23  (dimensions: 
60  m  length,  0.25  mm  ID,  0.25  um  film). 

Throughout  the  experimentation,  the  GC  was  operated  at  constant  conditions  (carrier:  helium  gas  with  flow  rate 
24  cm/sec;  oven:  140-240°C  at  4°C/min  and  injector:  250°C).  To  determine  glyceride  composition,  HPLC  (Agilent  1100 
Series  model)  equipped  with  refractive  index  detector  (model  G1362A),  Phenogel  100  A°,  300  x  7.8  mm  0.5  |jm  column 
(model  G1316A)  was  used.  The  calibration  of  both  the  analytical  instruments  has  been  done  with  authentic  standards. 

RESULTS  AND  DISCUSSIONS 
Characterization  of  Fish  Oil 

The  different  characteristics  of  tuna  fish  oil  were  found  as  density:  928  kg/m3;  acid  value:  0.507  mg  KOH/g  oil; 
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saponification  value:  180.0  mg  KOH/g  oil  and  iodine  value:  163.0  g  I2/100  g  oil.  The  analysis  with  HPLC  indicated  that 
tuna  fish  oil  contains  in  total  90  wt%  of  triglycerides.  According  to  the  GC  analysis,  the  tuna  fish  oil  was  found  to  contain 
mainly  fatty  acids  such  as  2.7  wt%  myristic  acid  (C14:0),  16.1  wt%  palmitic  acid  (C16:0),  3.6  wt%  palmitoleic  acid 
(C16:l),  4.03  wt%  stearic  acid  (C18:0),  13.1  wt%  oleic  acid  (C18:l),  1.3  wt%  linoleic  acid  (C18:l),  5.2  wt% 
eicosapentaenoic  acid  (C20:5),  26.0  wt%  docosahexaenoic  acid  (C22:6)  and  29.97  wt%  other  fatty  acids. 

Characteristics  of  Candida  rugosa  Lipase 

The  Candida  rugosa  lipase  used  has  got  the  following  characteristics.  Pore  diameter  is  48.7  (A°),  single  point 
total  pore  volume  is  0.09  (cc/g)  and  BET  surface  area  of  liapse  is  7.8  (m2/g)  [31].  The  higher  surface  area  of  lipase 
represents  the  better  chances  of  substrate  collision  to  the  binding  sites  present  on  the  surface  of  the  lipase.  Therefore 
increasing  the  number  of  collisions  per  second  increases  the  rate  of  the  reaction  [32].  The  shape  or  conformation  of  lipase 
active  site  plays  an  important  role  in  the  specificity  of  lipase  for  its  substrate. 

Therefore,  the  perfect  binding  of  substrate  depends  upon  the  uniformity  in  the  shape  and  size  of  the  lipases 
because  the  uniform  particles  of  lipase  could  only  provide  better  binding  of  substrate  in  the  open  lid  structure  of  lipase  in 
their  active  sites.  The  KM  constant  by  Michaelis-Menten  kinetics  is  a  mathematical  measure  to  represents  the  specificity  of 
lipases.  According  to  the  Paiva  et  al.  [32]  the  composition  of  the  lid  that  covers  the  active  site,  the  geometry  of  the  catalytic 
triad,  structure  and  the  dynamics  of  the  active  site  lid  opening  are  the  features  that  impart  lipases  their  unique  structure - 
function  characteristics. 

Effect  of  Reaction  Parameters 

The  effect  of  parameters  such  as  pH,  temperature,  speed  of  agitation,  water  concentration  and  solvent  amount  on 
the  activity  of  CRL  for  hydrolysis  of  tuna  fish  oil  in  biphasic  solvent  system  are  shown  in  Figures  1-5.  According  to 
Figure  1,  CRL  has  shown  maximum  activity  of  140.6  |imoles  FFAs  formed/ml  at  pH  6.5  while  the  activity  remains  129 
|imoles  FFAs  formed/ml  at  pH  7. 

It  shows  that  the  hydrolysis  reaction  can  be  carried  out  in  the  neutral  range  of  pH  6.5-7.0  because  lipases 
normally  contain  amino  acid  residue  as  an  active  site  and  it  generates  maximum  binding  affinity  for  the  substrate  in  the 
neutral  pH  range.  Therefore,  for  the  rest  of  the  study  it  has  been  selected  as  optimal  pH  6.5  for  CRL.  Han  and  Rhee  [33] 
have  reported  optimum  range  of  pH  value  i.e.  6.5  to  7.1  for  hydrolytic  reaction  of  olive  oil  catalyzed  by  CRL  enzyme. 
Similarly,  Knezevic  [28]  and  Okada  [34]  have  reported  pH  7.0  as  optimum  for  hydrolysis  of  palm  oil  and  sardine  oil 
catalyzed  by  CRL  enzyme. 

Maximum  activity  of  137.5  |imoles  FFAs  formed/ml  have  been  obtained  for  CRL  at  the  temperature  of  35°C  in 
1  h  as  shown  in  Figure  2.  CRL  lipase  is  vulnerable  at  the  higher  temperature  (>50°C)  indicating  a  sharp  reduction  in  the 
enzyme  activity  with  increase  in  the  temperature  because  it  results  in  the  denaturing  of  the  lipase  proteins.  Han  and  Rhee 
[33]  have  reported  optimum  temperature  in  the  range  of  30  to  35°C  for  hydrolysis  reaction  of  olive  oil  catalyzed  by  CRL 
enzyme.  Shimada  et  al.  [35]  have  used  35°C  temperature  for  hydrolysis  of  tuna  oil  by  CRL  enzyme.  Kahveci  et  al.  [36] 
have  also  reported  37°C  as  optimum  temperature  for  hydrolysis  of  salmon  oil  catalyzed  by  CRL  enzyme. 

The  effect  of  speed  of  agitation  was  studied  to  examine  the  mass  transfer  and  diffusion  resistances  on  the 
hydrolysis  reaction  of  oil  as  shown  in  Figure  3.  CRL  was  shown  to  give  maximum  product  formation  at  800  rpm 
equivalent  to  141  |amoles  FFAs  formed/ml  in  1  h.  It  is  clear  from  the  trend  observed  in  Figure  3,  that  initially  with  increase 
in  speed;  the  product  formation  was  increased  up  to  800  rpm.  Further  increase  in  the  speed  up  to  1000  rpm  did  not  lead  to 
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any  significant  improvement  in  the  product  formation,  indicating  that  reaction  is  proceeding  without  mass  transfer 
resistance  at  and  beyond  800  rpm.  According  to  Figure  4,  at  10:1  (w/w)  water  to  oil  ratio,  CRL  was  found  to  give 
maximum  activity  i.e.  138  umoles  FFAs  formed  per  ml  reaction  mixture  in  1  h.  Therefore,  hydrolysis  reaction  of  tuna  fish 
oil  with  CRL  was  carried  out  at  10:1  (w/w)  water  to  oil  ratio  to  obtain  best  results.  The  amount  of  water  used  plays  a 
crucial  role  in  hydrolysis  reaction.  In  case  of  oil  hydrolysis,  water  molecules  act  on  the  ester  bond  of  triglycerides  and 
provide  free  fatty  acids  and  glycerol  derivatives  as  product.  Similarly,  the  excess  amount  of  water  could  also  inhibit  the 
hydrolysis  reaction  by  competitive  inhibition  with  substrate  for  binding  to  enzyme.  CRL  enzyme  showed  maximum 
productivity  equivalent  to  136.6  pinoles  FFAs  formed  per  ml  in  1  h  at  1:1  (w/w)  solvent  to  oil  ratio  as  shown  in  Figure  5. 

It  shows  that  the  hydrolysis  activity  of  lipase  was  greatly  affected  by  the  concentration  of  iso -octane  as  solvent. 
Solvent  acts  as  an  organic  phase  in  the  biphasic  reaction  system  and  provides  the  ease  to  separate  product  from  the  reaction 
mixture.  Product  formed  after  the  reaction  has  tendency  to  be  soluble  in  solvent  and  therefore  are  easily  separable  due  to 
the  weak  interaction  with  enzymes  [37]  but  the  excess  amount  of  solvent  reduces  the  accessibility  of  enzyme  active  sites  to 
substrate.  Therefore,  the  most  suitable  reaction  conditions  for  hydrolysis  of  tuna  fish  oil  with  CRL  in  solvent  system  were 
obtained  as  pH  7.0,  temperature  35°C,  oil  to  water  ratio  of  1:10  (w/v)  and  oil  to  solvent  ratio  of  1:1  (w/v).  The  activity  of 
lipase  for  optimized  parameters  was  found  to  vary  from  137.5  to  141.2  for  CRL  in  pinoles  of  FFAs  formed  per  ml,  with 
results  reproducible  within  <+3%.  At  the  most  suitable  reaction  conditions,  86.5%  hydrolysis  was  achieved  with  CRL  in 
24  h.  Han  and  Rhee  [33]  have  also  reported  similar  suitable  reaction  conditions  for  the  hydrolysis  of  olive  oil  catalyzed  by 
CRL  enzyme  in  the  range  of  pH  6.5-7.1,  temperature  30-35°C  and  0.72-9.78  (v/v)  oil  to  water  ratio. 

Reaction  Mechanism 

The  lipase  catalyzed  hydrolysis  of  triglycerides  is  a  complex  reaction  that  involves  numerous  species  (mono-,  di-, 
and  triglycerides)  and  takes  place  by  consecutive  and  reversible  steps.  This  complexity  makes  it  difficult  to  obtain  a 
suitable  rate  equation  for  the  description  of  this  phenomenon.  Usually,  the  several  steps  (hydrolysis  of  mono-,  di-,  and 
triglycerides)  can  be  lumped  together  in  order  to  simplify  the  kinetic  analysis: 

Triglyceride  +  3H20  <-»  glycerol  +  3  fatty  acid  (1) 

Additionally,  inhibition  phenomena  by  fatty  acids  may  appear  as  the  reaction  proceeds.  The  enzyme  will  bind  to 
this  interface  changing  its  geometrical  conformation,  thus  exposing  the  active  center  and  making  it  accessible  to  the 
substrate  [39].  The  structure  of  the  active  center  of  lipases,  whose  high-resolution  three-dimensional  structures  have  been 
determined,  revealed  a  catalytic  triad  of  subcenters  Ser-His-Glu  for  Candida  rugosa  lipase,  which  are  involved  in  the 
binding  of  the  substrate  and  in  the  scission  of  the  ester  bonds.  The  existence  of  these  three  subcenters  may  also  justify  the 
appearance  of  higher  order  inhibition  in  the  case  where  the  inhibitor  attaches  itself  directly  onto  the  active  site  of  the 
activated  enzyme.  To  explain  the  results  obtained  in  this  work  and  taking  into  account  the  relatively  high  water  content,  a 
model  given  by  Prazeres  et  al.  [1],  including  nonlinear  product  inhibition  was  used,  considering  that  the  water  is  present  in 
excess  and  assuming  a  competitive  inhibition  of  lipase  by  the  fatty  acids  produced  upon  hydrolysis  of  the  triglycerides.  The 
lipolysis  mechanism  can  then  be  expressed  by  the  following  set  of  equations: 


KL  It, 

E  +  S1  <h>  ES1  ->  E  +  P 


(2) 


KB1 

E  +  P  «-»  EP 


(3) 


EP  4-  P  «->  EP2 


(4) 
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EP2+P^EP3  (5) 

The  number  of  product  molecules  that  we  assumed  could  be  bound  to  the  enzyme  was  arbitrarily  taken  as  3.  The 
value  of  3  corresponds  to  the  number  of  subsites  of  the  active  site,  derivation  of  the  rate  equation  from  the  steps  in 
equations  (2)  to  (5),  assuming  all  other  steps  to  be  in  equilibrium,  yields 

V  =   ^   (6) 

KMCl+Kj1P+Kj2P2+KjaPaO+S 

VM  =  K2Et  (7) 

where  V  =  initial  reaction  velocity  (pmoles  FFA/mlh),  VM  =  maximum  reaction  velocity  (pmoles  FFA/mlh), 
K2  =  enzymatic  rate  constant  (pmoles  FFA/mg  enzyme. h),  KM  =  Michaelis  Menten  constant  (pmoles  FFA/ml),  Kit  =Kel, 
Kj2  =  KeiKj.2  and  Ki3  =  KeiKe2Ke3 ,  Kn,  Ki2  ,  Ki3  =  first  order,  second  order  and  third  order  inhibition  constants  (p  moles 
FFA/mg  enzyme- h),  S  =  substrate  concentration  (pmoles  FFA/ml),  P  =  product  concentration  (pmoles  FFA/ml),  E,=  total 
enzyme  concentration  (mg/ml). 

Application  of  Prazeres  Model  for  Hydrolysis  Kinetics  of  Tuna  Fish  Oil  Using  CRL 

The  model  equation  given  by  Prazeres  et  al.  [1]  as  shown  in  Equation  6  was  used  in  the  present  work  considering 
second  order  product  inhibition  for  which  the  rate  equation  is  as  follows: 

dP  K?EtS 
V=  —  = 


dt     ^(l+K^P  +  K^P^  +  S  (8) 

The  above  model  was  used  to  fit  the  experimental  data  to  determine  the  kinetic  constants  and  also  the  values  of 
product  (FFAs)  formed  with  reaction  time  by  conducting  MATLAB®  simulation  for  four  different  initial  substrate 
concentrations  independently.  The  comparison  of  predicted  and  experimental  data  of  product  formation  with  respect  to 
time  of  reaction  is  shown  in  Figure  6.  The  kinetic  constants  obtained  by  fitting  model  for  different  initial  substrate 
concentrations  independently  using  MATLAB®  are  given  in  Table  1.  The  order  of  the  values  of  KM  (2.31  to  6.62  pmoles 
FFA/ml),  VM  (297.5  to  675.2  pinoles  FFA/ml -h),  Ku  (2.5  X  10"6  to  9.5  X  10"7  pmoles  FFA/mg  enzyme-h),  Kl2  (0.034  to 
0.06  pinoles  FFA/mg  enzyme-h)  and  K2  (74.4  to  168.8  pinoles  FFA/mg  enzyme-h)  obtained  by  using  Prazeres  model  is 
almost  similar  for  different  initial  substrate  concentrations. 

The  comparison  of  experimental  and  predicted  values  of  product  formed  with  respect  to  reaction  time  obtained  by 
MATLAB®  simulation  taking  all  the  four  sets  of  data  together,  are  shown  in  Figure  7  and  the  value  of  kinetic  constants 
were  obtained  as  shown  in  Table  1.  The  values  of  R2  has  been  found  in  the  range  of  0.952  to  0.973  and  the  values  of 
RMSE  has  been  found  in  the  range  of  43.05  to  126.7,  indicating  a  close  fitting  between  the  predicted  and  the  experimental 
data. 

CONCLUSIONS 

The  most  suitable  reaction  conditions  for  hydrolysis  of  tuna  oil  with  CRL  in  biphasic  solvent  system  have  been 
determined.  The  lipase  was  observed  to  work  optimally  at  pH  7,  temperature  35 °C,  l:l(w/w)  oil  to  solvent  ratio  and  1:10 
(w/w)  oil  to  water  ratio.  For  hydrolysis,  the  kinetic  model  given  by  Prazeres  et  al.  [1],  incorporating  second  order  product 
inhibition  was  used.  The  values  of  kinetic  constants  have  been  evaluated.  Two  statistical  parameters  R2  and  RMSE  were 
estimated  to  check  the  goodness  of  fitting  between  the  predicted  and  experimental  values  of  product  formation  with  respect 
to  time  of  reaction. 
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NOMENCLATURE 

CRL=  Candida  rugosa  lipase 

DDGS  =  Distillers  dried  grain  with  solubles 

DHA  =  Docosahexaenoic  acid 

FFA  =  Free  Fatty  acid 

Et  =  Total  enzyme  concentration  (mg/ml) 

K2  =  Enzymatic  rate  constant  (|i  moles  FFA/mg  enzyme  h) 

KM  =  Michaelis  Menten  constant  (\x  moles  FFA/ml) 

Ku,  Ki2,  Ki3  =  first  order,  second  order  and  third  order  inhibition  constants  (|imoles  FFA/mg  enzyme  h) 

P  =  Product  concentration  (|i  moles  FFA/ml) 

PUFA  =  Poly  Unsaturated  Fatty  acid 

S  =  Substrate  concentration  (p,  moles  FFA/ml) 

t  =  time  (h) 

V  =  initial  reaction  velocity  (umoles  FFA/ml  h) 
VM  =  maximum  reaction  velocity  (umoles  FFA/mlh) 
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APPENDICES 


Table  1:  Kinetic  Constant  at  Individual  and  Combined  Substrate  Concentrations 
during  Hydrolysis  of  Tuna  Oil  with  CRL  Using  Prazeres  Model 


Substrate 
Concentrations 
(umoles/ml) 

KM 
(umoles 
FFA/ml) 

VM 
(umoles 
FFA/ml- h) 

K„ 

(fimoles  FFA/mg 
enzyme- h) 

Ki2 
(umoles 
FFA/mg 
enzyme- h) 

K2 
(umoles 
FFA/mg 
enzyme- h) 

R2 

RMSE 

S01  214 

2.31 

297.5 

9.5  X  10"7 

0.060 

74.4 

0.961 

43.0 

S02  267 

2.44 

318.0 

2.5  X  10"6 

0.039 

79.5 

0.973 

45.6 

S03  334 

5.69 

666.8 

1.0  X  10"' 

0.036 

166.7 

0.952 

68.5 

S04401 
Average 

6.62 
4.26 

675.2 
489.4 

3.3  X  10"6 
6.0  X  10"6 

0.034 
0.042 

168.8 
122.3 

0.960 
0.961 

69.5 
56.7 

Using  all  the 
four  sets  of  data 
together 

9.90 

503.6 

6.0  X  10"4 

0.017 

125.9 

0.957 

126.7 
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Figure  1:  Effect  of  pH  on  Activity  of  CRL  during  Hydrolysis  of  Tuna  Fish  Oil 
{Experimental  Conditions:  1  g  Tuna  Oil,  1  ml  Iso-Octane,  10  ml  Water,  Temperature  35°C, 
Agitation  Speed  800  rpm,  40  mg  Enzyme,  Reaction  Time  1  h} 
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Figure  2:  Effect  of  Temperature  on  Activity  of  CRL  during  Hydrolysis  of  Tuna  Fish  Oil 
{Experimental  Conditions:  1  g  Tuna  Oil,  1  ml  Iso-Octane,  10  ml  Water,  pH  6.5, 
Agitation  Speed  800  rpm,  40  mg  Enzyme,  Reaction  Time  1  h} 
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Figure  3:  Effect  of  Agitation  Speed  (rpm)  on  Activity  of  CRL  during  Hydrolysis  of  Tuna  Fish  Oil 
{Experimental  Conditions:  1  g  Tuna  Oil,  1  ml  Iso-Octane,  10  ml  Water,  pH  6.5, 
Temperature  35°C,  40  mg  Enzyme,  Reaction  Time  1  h} 
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Figure  4:  Effect  of  Oil  to  Water  Ratio  (w/v)  on  Activity  of  CRL  during  Hydrolysis  of  Tuna  Fish  Oil 
{Experimental  Conditions:  1  g  Tuna  Oil,  1  ml  Iso-Octane,  pH  6.5,  Temperature  35°C, 
Agitation  Speed  800  rpm,  40  mg  Enzyme,  Reaction  Time  1  h} 
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Figure  5:  Effect  of  Oil  to  Solvent  Ratio  (w/v)  on  Activity  of  CRL  during  Hydrolysis  of  Tuna  Fish  Oil 
{Experimental  Conditions:  1  g  Tuna  Oil,  10  ml  Water,  pH  6.5,  Temperature  35°C, 
Agitation  Speed  800  rpm,  40  mg  Enzyme,  Reaction  Time  1  h} 
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Figure  6:  Comparison  of  Experimental  (1— >0,  2— >A,  3— >o,  4— >+)  and  Predicted  (1- 
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 )  Data  of  Product  Formation  at  Initial  Substrate  Concentration  SOI,  S02,  S03  and  S04  Using  Prazeres  Model 

{Experimental  Conditions:  1  g  Tuna  Oil,  10  ml  Water,  Isooctane  in  Suitable  Amount,  pH  6.5,  Temperature  35°C, 
Agitation  Speed  800  rpm,  40  mg  Enzyme,  Batch  Volume  12  ml} 
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Figure  7:  Comparison  of  Experimental  (1— >0,  2— >A,  3— *o,  4— >+)  and  Predicted  (1- 
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 )  Data  of  Product  Formation  after  Combining  Four  Data  sets  S0i,  S02,  S03  and  S04  Together  Using  Prazeres 

Model  {Experimental  Conditions:  1  g  Tuna  Oil,  10  ml  Water,  Isooctane  in  Suitable  Amount,  pH  6.5, 


Temperature  35°C,  Agitation  Speed  800  rpm,  40  mg  Enzyme,  Batch  Volume  12  ml} 


